Gate control of spin dynamics in III— V semiconductor quantum dots 
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We show that the g-factor and the spin-flip time T\ of a heterojunction quantum dot is very 
sensitive to the band-bending interface electric field even in the absence of wave function penetration 
into the barrier. When this electric field is of the order of 10 V/cm, g and Ti show high sensitivity 
to dot radius and magnetic field arising from the interplay between Rashba and Dresselhaus spin- 
orbit interactions. This result opens new possibilities for the design of a quantum dot spin quantum 
computer where g-factor and Ti can be engineered by manipulating the spin-orbit coupling through 
external gates. 

PACS numbers: 73.21.La; 03.67.Lx; 71.70.Ej; 85.35.Be; 
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Understanding and controlling the behavior of spins 
in semiconductor heterostructures may lead to a whole 
new class of devices^* 2 " ranging from spin-polarized p-n 
junctions 3 " to a quantum dot spin quantum computer, 4 
Electrical control over spin-orbit coupling parameters in 
a quantum well has long been suggested as an effec- 
tive means to manipulate spin^ and was recently exper- 
imentally demonstrated in InAs heterostructures^ This 
was possible because the Rashba^ and the Dresselhaus 9 
spin-orbit interactions are sensitive to the electric field 
providing vertical confinement to a two dimensional elec- 
tron gas (2DEG), which is approximately proportional 
to the 2DEG density and can be controlled through a 
gate voltage. When additional gates provide lateral con- 
finement, a single electron can be trapped in a quan- 
tum dot (QD), 10 whose orbital states in an external 
magnetic field perpendicular to the 2DEG are the well 
known Fock-Darwin statesiii Here we consider the effect 
of Rashba and Dresselhaus spin-orbit interactions in the 
spin-doublet ground state of a Fock-Darwin QD. Using an 
effective mass approximation and exact diagonalization 
of a Fock-Darwin subspace we show that the ground state 
g-factor and spontaneous phonon emission rate \jT\ 
(due to the spin-orbit admixture mechanism) 1 ^ can be 
substantially manipulated by varying the heterojunction 
electric field in the range of 10 5 — 10 6 V/cm (correspond- 
ing to 2DEG density of the order of 10 12 - 10 13 cur 2 ). 

Recently, electrical control over GaAs quantum well 
g-factor has been achieved by forcing the electron wave 
function to overlap with the AlGaAs barrieriA&ii Here 
we intentionally neglect barrier penetration, to show that 
overlap with a different material is not a necessary con- 
dition to achieve substantial electrical control over quan- 
tum dot g-factor. Moreover, by avoiding barrier pen- 
etration one can suppress an additional spin-lattice re- 
laxation mechanism due to interface motioni^ These re- 
sults open new possibilities in the design of a quantum 
dot quantum computer: for example, spin qubits can be 
brought in and out of resonance to a global spin reso- 
nance field by gate control of their g-factor;— a speed up 
in quantum computer initialization (setting all spins up 
for example) can be achieved by decreasing the spin-flip 



time T\ with a gate voltage. 

The Hamiltonian for a single electron bound to an het- 
erojunction quantum dot can be divided into five parts, 



H = Ho + H z + Hr + "Hx>i + Ht>2- 



(1) 



The first contribution corresponds to a single 2D electron 
confined in the xy plane by a parabolic potential and 
subject to a magnetic field B, 



Ho 



2m* 



(2) 



where the kinetic momentum P = p + e/cA is written 
with the canonical momentum p = —ih(d x ,d y ,0) and 
vector potential A = B/2(—y,x,0) confined to the 2D 
plane. Here e is the electron charge, c is the velocity of 
light, m* is the conduction band edge effective mass, wo is 
the parabolic confining potential frequency, go is the Bulk 
g-factor, fig is the Bohr magneton, and o~ z is the diagonal 
Pauli matrix. Ho is diagonal when written as a function 
of the Fock-Darwin number operators n± = a!j_a±^i 



Ho = Huj+ n 
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-( X ±iy)--(d x ±id y ), (4) 

1 I 

— (it iy) + - (d x Tid y )- (5) 



Here u>± = f2 ± ui c /2, with = a/Wq + u c /4 an d 
lu c = eB/m*c being the renormalizcd dot frequency and 
cyclotron frequency respectively, with £ — y/h/m*£l be- 
ing the Fock-Darwin radius which sets the length scale 
for the eigenstates \n + n_a) (a = ±1 represents the spin 
up/down states in the z direction), The second term in 
the Hamiltonian Q represents the quantum well confine- 
ment in the growth (z) direction, H z — p 2 z /2m* + V(z), 
where V(z) is a triangular well, V(z) — eEz for z > 
and V(z) — oo for z < 0. A simple numerical calculation 
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FIG. 1: Quantum dot g-factor divided by the bulk g-factor 
as a function of the dot radius lo = v/7i/m*a>o f° r a GaAs 
heterojunction. Each curve corresponds to a different electric 
field: From top to bottom, E = 10 4 , 1, 2, . . . , 10 x 10 5 . The 
magnetic field applied in the growth [100] direction is assumed 
to be 1 Tesla. At large lo the two lowest energy states have 
same spin, hence the level crossing. The results for GaSb are 
similar, except that the level crossing occurs for smaller lo- 

leads to the H z ground state 

*o.(*) = 1.4261k 1 / 2 Ai(/w + Ci), (6) 

where £i = —2.3381 is the first zero of the Airy function 
Ai, while the inverse length scale k is set by 

k = (2m*eE/h 2 ) 1/3 , (7) 

and the ground state energy is E Qz — —Q±eE/K. In the 
discussion below we will make use of the average mo- 
mentum squared in the state |JB}, (p 2 ) = 0.7794(fi,K) 2 , 
and the average position (z) = 1.5587/k (which is the 
thickness of the 2DEG). We now turn to the spin-orbit 
interactions, third to fifth terms in Eq. QJ. A k • p band 
structure calculation for zincblende materials^ leads to 
the bulk conduction band spin-orbit interaction 

H B uik = Jc/(2h 3 )a ■ P, (8) 

where P x — P X (P 2 — P 2 )+h.c, P y and P z can be obtained 
by cyclic permutations. Notice that Eq. JSJ is hermitian 
and gauge invariant. The value of j c is determined by the 
band structure parameters of the III-V semiconductors 
(Table I). By averaging Eq. JSJ) over the quantum well 
ground state [Eq. ©] we get two spin-orbit terms, linear 
and cubic in momenta (here the quantum well growth 
direction is assumed to be [001])ji£ 

Hoi = Q.77M lc K 2 /h{-a x P x +<r y P y ) : (9) 

H D2 = lc/h 3 (a x P X P 2 - GyPyPl) + h.C. (10) 

The structural inversion asymmetry in V{z) leads to the 
Rashba interaction^ 

Hn = a R eE/h{u x Py-a y P x ), (11) 
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FIG. 2: Quantum dot g-factor as a function of the magnetic 
field B for a GaAs heterojunction. The electric fields are the 
same as in Fig. 1. Dot radius is lo = 20 nm. At large B a 
similar level crossing as in Fig. 1 takes place. 

which for the triangular well considered here is di- 
rectly proportional to E (a R depends on band structure 
parameters) £ It is useful to condense Hn and Hui m 
a single Hamiltonian written as a function of the Fock- 
Darwin operators, 

Hbi+H r = V(t + + V^<t-, (12) 



V = — ct'_aL + Q'+a + + i/3_a_ — ift + cJ + , (13) 

with a± = (a x ± ia y )/2 and spin-orbit energy scales de- 
fined as 

ot± = a R eE£,±, (14) 
f3 ± = 0.78 7cK 2 C±, (15) 
1 eBt 

£t = 7 ± 2fe' (16) 

Therefore Hr and Hui couples Fock-Darwin levels dif- 
fering by one quantum number. The amount of spin 
up/down admixture in the ground state is given by the 
ratios a±/hwQ, 0±/hoJo, which at B w are directly 
proportional to the dot radius l = \Jh/m*ujQ. There- 
fore larger dots will be more sensitive to spin-orbit cou- 
pling, at least within perturbation theory (see below). 
Moreover, the question of whether Rashba or Dresselhaus 
dominates depends only on material parameters and elec- 
tric field, since a±/(3± = a R / lc {h 2 /m*) 2 / 3 {eE) 1 / 3 . If 
the electric field ranges from 10 4 — 10 6 V/cm this ratio 
equals 0.1 -0.7 for GaAs, 0.2 -1.0 for GaSb, 1.5-6.8 for 
InAs, and 5.6 — 26 for InSb. Therefore for III-V semicon- 
ductor quantum dots it is important to consider the in- 
terplay between Rashba and Dresselhaus spin-orbit inter- 
actions. High electric field GaSb heterojunctions might 
realize the condition a± = (3±, leading to an interesting 
simplification of Eq. i|12|) , which becomes proportional to 
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(a x + a y )ML However this symmetry is broken by Eq. J3J, 
which contains a magnetic field pointing in the [001] di- 
rection. Hence it will have no consequences here [the 
highly symmetric case of B | [110] and a± = j3± leads 
to QD ^-factor exactly equal to the bulk value go, with 
no spin-lattice relaxation, as long as the cubic spin-orbit 
interaction Eq. I|10|) is neglected]. 

We now turn to the cubic spin-orbit term [Eq. lfTU|) ]. 



\\a} 2 a + + A2fl_ l | Z a_ + A4<2_ 



t 2 , 



-A 3 o^(n + + l) + ^a_(2n + + l) 



3 +v 

,t3 



\4c1_ + Asa^_ 



1 \ , ^2 +2 t ^1 2 

1) + y a + a - ~ y a + a - 

+h.c, 



(17) 



with Ax - 3/4 7c £+e, A 2 = 3/4 7c C-^, A 3 = l/4rf£%, 
A 4 = l/4 7c £?.. Hr>2 is often neglectedj±£ a well justified 
approximation for heterojunctions with a small Fermi 
wave vector (H.D2 was considered recently in a different 
context) However, we will show that 7Yd2 leads to two 
interesting effects for small few electron quantum dots. 
The first term in Eq. I|17|) with its hermitian conjugate 
couples the state \n + n^a) with \n + + a, n_ — 2a, — a) 
which are degenerate at lu c « (the exact location 

of the anticrossing depends on Zceman splitting). The 
magnitude of the anticrossing is given by 



AE: 



V3 ; 3 



v/n+(n_ + l)(n_ + 2). 



(18) 



For a GaAs dot with l = 10 nm, AE ~ 0.1 meV, while 
for GaSb, In As, and InSb it can reach 1 meV. Note, how- 
ever, that this anticrossing appears only for n + +rt_ > 1. 
It may have interesting consequences for spin dependent 
transport through a few electron QD. The third to sixth 
terms in Eq. I|17|) are linear in a±, a± leading to an en- 
hancement of the spin-orbit effect [Eq. (|13fl ] . This affects 
our g-factor and Xf calculations by as much as a factor 
of 2, when Iq < 10 nm. 

Our QD ground state g-factor is defined by g = (E 2 — 
Ei)/(hbB), Ei and E 2 being the ground and first excited 
states including spin. Considering Eq. i|12|) as a second 
order perturbation to Eq. with B field in the [001] 
direction we get 



9o 



h A 

1 m„m* 3 



[0.6 7c V(l - §) - a 2 R e 2 E 2 (l + 5)] I 



{0.6 7 2 c k 4 (1 - 5 + 8 2 + S 3 ) 



2 h 6 

-a 2 R e 2 E 2 (l + 5 + S 2 -S 3 )}u; 2 Jt 

+o ( Wc M>) 4 . 



(19) 



Here 5 = gom*/m e , this expression being valid up to 
second order in lu c /ujq and the spin-orbit admixtures. 
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FIG. 3: Ratio between dot g-factor and bulk (?-factor as a 
function of the dot radius Iq for an InAs heterojunction. Each 
curve corresponds to an electric field, from bottom to top, 
E = 10 4 , 1, 2, . . . , 10 x 10 5 . The magnetic field applied in 
the growth [100] direction is assumed to be 1 Tesla. Notice 
the qualitative difference with respect to Fig. 1: Here Rashba 
interaction dominates the g-factor, while in Fig. 1 Dresselhaus 
dominates. The results for InSb are similar. 
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FIG. 4: j/jo as a function of magnetic field for an InAs het- 
erojunction quantum dot with io = 20 nm. Electric fields are 
the same as in Fig. 3. The results for InSb are similar. 



Clearly one sees that <?-factor displays a rich behavior 
as a function of QD radius and electric field. In particu- 
lar, g — go will be positive and proportional to E 4 ^ 3 if the 
Dresselhaus spin-orbit interaction is dominating (GaAs 
and GaSb), but negative and proportional to E 2 when 
Rashba dominates (InAs, InSb). In addition these ef- 
fects increase with increasing dot radius lo, and there is 
a B 2 dependence at higher magnetic fields. 

We study the g-factor behavior at high electric and 
magnetic fields, and large dot radius by resorting to 
exact diagonalization of the full Hamiltonian [Eq. QJ] 
(similar calculations for Landau levels in a quantum well 
are available). 20 Our basis consists of nine Fock-Darwin 



shells (n_| 



< 



), which together with spin leads 
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to a 90 x 90 matrix for the Hamiltonian. Figs. 1 and 
2 show results for GaAs, Figs. 3 and 4 for InAs. GaSb 
shows very similar behavior to GaAs, while InSb is sim- 
ilar to InAs. The differences between the two sets of 
materials is attributed to Dresselhaus spin-orbit interac- 
tion dominating in Figs. 1 and 2, but Rashba dominating 
in Figs. 3 and 4 (the Rashba interaction appears to be 
dominant in SiGe heterostructures, therefore g-factor be- 
havior should be similar to the InAs case considered here, 
except that ckr is three orders of magnitude smaller and 
g » 2; hence the effects considered here should be quite 
small for Si heterojunctions: we estimate g — go ~ — 10~ 3 
for E ~ 10 5 V/cm)j2i Our result suggests g-factor can 
be controlled by a gate voltage (that either changes the 
longitudinal electric field E or the lateral confinement 
Iq) as long as E ~ 10 5 V/cm. An important feature 
of Figs. 1-4 is that a level crossing takes place for large 
enough Iq and B. In this regime the two lowest energy 
states of the QD are approximate Landau levels with 
the same spin, leading to extremely fast phonon emis- 
sion rates (see below). Hence a QD quantum computer 
should operate away from this level crossing, which is ac- 
tually a smooth anti-crossing for InAs (see Figs. 3 and 
4). Note that Figs. 1-4 plot the ratio of QD to bulk 
g-factor (g/go), therefore the corresponding deviation in 
QD Zceman energy from the bulk value will be rather ap- 
preciable for GaSb, InAs, and InSb since these materials 
have g 10 (Table I). 

We now turn to calculations of the transition rate 
between the two lowest energy states due to spon- 
taneous phonon emission. The electron-piezophonon 
interaction 12 



U" 



- P h 



pi(q-r-w, Q t) A Jkt 

2pVcu qa e 6 qa ^ 



h.c. 



(20) 
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FIG. 5: Spin-lattice relaxation rate for the spin-doublet 
ground state of a 20 nm one-electron quantum dot as a func- 
tion of the external magnetic field. Notice for all materials 
1/Ti oc B 5 at low B fields. At high B the Zeeman phonon 
wavelength becomes smaller than the dot radius and the rate 
is strongly suppressed. This effect is evident at B ~ 1 Tesla 
for the narrow gap materials, which have quite large bulk 
o-factors. Here we assumed E = 10 4 V/cm. The curves 
terminate at the level crossing, when the rate is extremely 
enhanced because the transition states have the same spin. 
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GaAs GaSb InAs InSb 
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m* /m e 



-0.44 -7.8 -15 -50.6 
0.067 0.0412 0.0239 0.0136 
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[10 5 cm/s] 


3.03 
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2.35 


2.29 


P 


[g/cm 3 ] 


5.3176 


5.6137 


5.6670 


5.7747 



couples these states in the presence of spin-orbit admix- 
ture. Here b qa creates an acoustic phonon with wave 
vector q and polarization e Q (a = L, Ti,T%), p is the ma- 
terial density, and V the volume of the sample. A qa is 
the amplitude of the electric field created by the phonon 
strain, which is given by qiqke(3ijke qa , with q = q/q, 
sPijk = ehu (see Table I) for i ^ k, i ^ j, and j =/= k. 
The polarization directions are 

?l = (sin # cos </>, sin sin ^, cos 9), (21) 
en = (cos 9 cos <j>, cos 9 sin <j>, — sin 9) , (22) 
&T2 = (— sin<^>, co$,(f>, 0). (23) 

The transition rate is given by Fermi's golden rule, 



V 



Tl • ,X J d 3 q\{l\U q e % h \2)\ 2 S(hLJ qa -E 2 + E 1 ), 

. (24) 

which under the same perturbative approximation as 



TABLE I: Parameters used in our calculations. 8 " 



Eq. O leads to 
4 



1 



3 1 \ Am* 4 (eh M ) 2 



105tt V 4 4 4 



ph 7 



.2..] . 2 ,2ir2\ / \ 



x 0.6l7> 4 + ate'E' 



x(n ph + l) l + O(u c /iu f 



(25) 



Here st and sl are the transverse and longitudinal acous- 
tic phonon velocities respectively. The spin-flip rate is 
extremely sensitive to QD radius and external magnetic 
field. At temperatures lower than Zeeman splitting, the 
emitted phonon occupation number n p h is much smaller 
than one, and Eq. l (25Tl is independent of temperature. At 
higher temperatures Raman processes will dominated 

It is interesting to study deviations from the perturba- 
tive approximation Eq. I|25|) . In particular, at high mag- 
netic fields the resonant phonon wavelength Xz = hs/ 'Ez 
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becomes much smaller than the dot radius making the 
dipolar approximation on the electron-phonon interac- 
tion inappropriate [Eq. Q25JI assumes the exponent in 
Eq. 1)21) [) can be approximated by ~ 1 + iq • r]. Fur- 
thermore, one immediately sees that 1/Ti is extremely 
sensitive to the energy difference E1 — E2, assumed equal 
to the bulk Zeeman energy in Eq. I|25[l . Here we show cal- 
culations of Eq. I|24|) using energy levels and eigenstates 
obtained by exact diagonalization in a 90 dimensional 
Fock-Darwin basis. In addition, we go beyond the dipo- 
lar approximation by using the identity 



o*q r — 



(26) 



where 77(6*, <f>) — qis'mOe 1 ^ /2 depends on the polar angles 
of the phonon wave vector q [because of this dependency, 
we have to perform the angular integrals in Eq. (|24|l nu- 
merically] . Each of the exponents in Eq. I|26[) is expanded 
in powers of 77, but we note that within our subspace 
n ± < 8, therefore only up to the ninth power needs to be 
retained. We checked the convergence of our calculations 
by reducing the Fock-Darwin subspace and noting that 
no appreciable change takes place for n + + n_ > 4. Our 
results agree with perturbation theory [Eq. J2HJ)] at low 
B and E. Fig. 5 shows the spin-flip rate as a function 
of the magnetic field. It is evident that materials such 
as InAs and InSb deviate from perturbation theory by 
more than three orders of magnitude when B is as low as 
1 Tesla. This happens because taking into account the 
full electron-phonon Hamiltonian leads to an exponential 
decrease in the rate when ql ^> 1, since Eq. I|2(J|I oscil- 
lates appreciably in this regime. Note, however, that 
1/Ti cx B 5 at low enough B for all materials. Fig. 6 
shows the dependency of the spin-flip rate with lateral 
confinement radius Iq in a GaAs QD. As the electric field 
increases, the dependency with Iq displays an striking be- 
havior, which happens due to the sign change of <?-factor 
shown in Fig. 1. A small Zeeman energy implies neg- 
ligible phonon density of states, and hence the rate is 
zero for l ~ 50 nm and E = 7 x 10 5 V/cm in Fig. 6. 
We expect 1/Ti will behave similarly to Fig. 6 when g 
changes sign due to barrier penetration in AlGaAs^ this 
property being extremely useful in the initialization and 
decoherence suppression of a QD quantum computer. 

We now discuss possible corrections to the simple 
model discussed here. At strong confinement in the 100 
direction one expects T-X valley mixing to become im- 
portant. For GaAs, Erx — 0.48 eV, which is comparable 
to E 0z [Eq. ©] only when the electric field is the highest 
considered here, E > 10 6 V/cm. This also holds true for 
InAs and InSb, but in GaSb T-X coupling will be impor- 
tant for E > 5 x 10 5 V/cm. Therefore even though a full 
kp calculation would yield some corrections^ we do not 
expect it to change our results qualitatively in the range 
considered here. The same holds true for the inadequacy 
of the Rashba Hamiltonian, which starts deviating from 
Eq. JTTJ when E ~ 10 6 V/cm£ 

Before concluding, we discuss the approximations and 
the limitations of this work. The most essential approx- 




«>/„ 



1X10 V/cm 



5X10- 



7X10" 



150 

Quantum dot radius I. [nm] 



200 



FIG. 6: Spin-flip rate 1/Ti due to spin-orbit admixture as a 
function of dot radius in GaAs. For small radius, 1/Ti oc Iq. 
For E = 7 x 10 5 V/cm there is a striking change in behavior: 
This happens due to the sign change in g-factor seen in Fig. 1. 



imation of our model, the use of the k • p perturbation 
theory within an effective mass approximation scheme to 
describe the conduction band, has been extensively used 
in the literature 2 ^ and should be well-valid for the prob- 
lem we study. We have made two additional non-essential 
approximations in our theory in order to simplify our 
numerical computations: The triangular well approxi- 
mation for the z-confinement of the wave function and 
the parabolic well Fock-Darwin confinement approxima- 
tion in the 2D x — y plane. These approximations are 
reasonable enabling us to produce numerical results for 
a range of system parameters in several different semi- 
conductor structures, which would have been difficult, if 
not impossible, to carry out had we used more realis- 
tic (and thereby numerically more demanding) quantum 
dot confinement potentials (The triangular well approx- 
imation was employed recently to derive new results re- 
garding D'yakonov-PereP relaxation anisotropy for con- 
duction electron spins confined by heterojunctions with 
E > 10 5 V/cm)i2i Our most important qualitative result, 
establishing the viability of controlling the spin dynam- 
ics of III-V semiconductor quantum dots (both (/-factor 
and T\ engineering) by using external gates to suitably 
manipulate the spin-orbit coupling through Drcssclhaus 
and Rashba effects, should be completely independent 
of these approximations. In fact, we expect that the 
use of more sophisticated confinement models may ac- 
tually make the gate control effects we predict somewhat 
stronger by pushing the required electric fields (and con- 
sequently 2D carrier densities) to somewhat lower values 
than our predicted 10 5 V/cm range. The main limita- 
tion of our predicted spin-orbit coupling induced gate 
control effect is, in fact, the rather large electric fields 
(~ 10 5 V/cm) and the associated 2D carrier densities 
(~ 10 12 cm -2 ) that are required to produce significant 
gate control effects. 
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In conclusion, we show that quantum dot longitu- 
dinal g-factor and spin-flip time T\ can be controlled 
electrically even in the absence of wave function over- 
lap with a different material. These parameters show 
a striking dependence with dot radius and magnetic 
field when the 2DEG confinement is strong (electric field 
E ~ 10 5 V/cm). For example, the g-factor changes sign 
and T\ is extremely sensitive to the dot radius, ^-factors 
for one-electron dots can be measured using transport 
spectroscopy^ We show that Ti is drastically increased 
in narrow gap materials (InAs, InSb) due to deviations 
from the dipolar approximation in the electron-phonon 
interaction, suggesting these materials are promising for 
the fabrication of a quantum dot spin quantum computer. 
Ti is found higher than 10 -4 s under quite different cir- 
cumstances (see Figs. 5 and 6) showing that small III-V 
one-electron quantum dots (^o < 50 nm) will have their 
low temperature phase coherence time (T2) dominated 
by nuclear induced spectral diffusion. 26 This result es- 
tablishes the versatility of III-V quantum dots as units 
for single spin manipulation. A related finding of interest 



in our work is the dual importance of both Dresselhaus 
(i.e. the bulk inversion asymmetry inherent in Zincblcnde 
structures of III-V semiconductors) and Rashba (i.e. the 
real space structural inversion asymmetry present in a 
heterostructure due to external electric fields) spin-orbit 
coupling terms in semiconductor nanostructures - in par- 
ticular, for GaAs and GaSb quantum dot structures in- 
vestigated in this work, we typically find the bulk in- 
version asymmetry (i.e. Dresselhaus) effect to be quan- 
titatively more important than the Rashba effect. The 
relative quantitative importance of the Dresselhaus effect 
in III-V nanostructures should have considerable signif- 
icance not only in the g-factor engineering and the spin 
relaxation time control of relevance to the spin quantum 
computer architecture (that we consider in this work), 
but also in the fabrication of the Datta-Das spintronic 
transistoi^ where spin-orbit coupling is used to modu- 
late a spin-polarized current in a field effect transistor 
configuration. The authors acknowledge discussions with 
A. Kaminski and I. Zutic. This work is supported by 
ARDA, LPS, US-ONR, and NSF. 
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